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ABSTRACT 

A Landsat Thematic Mapper quality evaluation study was conducted to 
determine the existence of geometric and radiometric sensor errors in 
data in the post-launch environment. (The Landsat-4 Thematic Mapper 
will be denoted TM-4 and similarly for Landsat-5.) The study began with 
the launch of Landsat-4; several error conditions were found, including 
band-to-band geometric misregistration and detector-to-detector radiome- 
tric calibration errors. These results were reported (Anuta, 1984) and 
other investigators and used by NASA to update the system. 

This paper describes similar analysis for the Landsat-5 Thematic 
Mapper and comparisons are made with results for Landsat-4. Remaining 
band-to-band misregistration was found to be within tolerances and 
detector-to-detector calibration errors were not severe but still exis- 
tent. More coherent noise signals were observed in TM-5 than in TM-4, 
although the amplitude was generally less. The scan-direction differ- 
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ences observed in TM-4 were still evident in TM-5, the largest effect 
being in Band 4 where nearly a one digital count difference wa3 
observed. Resolution estimation was carried out using roads in TM-5 for 
the primary focal plane bands rather than field edges as in TM-4. Esti- 
mates using roads gave better resolution than the ones using edges. 
Thermal IR band calibration studies were conducted and new nonlinear 
calibration procedures were defined for TM-5. The overall conclusion is 
that there are no first-order errors in TM-5 and any remaining problems 
are second or third order. 


Introduction 

This paper describes analysis of Landsat 5 Thematic Mapper data to 
evaluate geometric and radiometric qual.ty. In a previous paper 

A 

(Anuta, 1984), Landsat 4 Thematic Mapper and multispectral scanner data 
were analyzed for the geometric and radiometric quality and Information 
content. The results of the Landsat 5 work are compared to that for the 
Landsat 4, so some of the material in the previous paper is repeated 
here. Not all the topics studied in the Landsat 4 study were repeated, 
so only a subset of the results are comparable. As a guide to the 
results in the two papers, the subjects are listed below, with those 
reported in this paper indicated in the second column. 


a 


See Appendix for Landsat-4 resulrs 
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IEEE 

GRS ASP 
May 8*) Sept. 

1. Effects of resampling (comparison of A & P tapes) X 

2. Band-to-band registration X X 

3. Detector- to-detector calibration X X 

4. Coherent-noise analysis.... X X 

5. Spatial resolution estimation X X 

6. Scan-angle effects... X 

7. Scan-direction effects X X 

8. Spectral dimensionality X 

9. Spectral separability X 

1 0 . Thermal calibration X X 


The results of the Landsat-4 study revealed several geometric and 
radiometric errors which, when reported to NASA by the authors as well 
as by several other investigators, resulted in corrections and modifica- 
tions that were incorporated in Landsat 4 ground processing and the 
Landsat 5 sensors. This study is a final evaluation of the success of 
these changes and the results indicate no significant problems with the 
TM system. Some second-order problems with noise and calibration exist 


and these are discussed. 
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Band-to-Band Misregistration 

The importance of band-to-band registration is that it specifies 
the spatial alignment of all spectral bands in the system. Exact 
requirements for band-to-band misregistration were as follows [ 3 ]: "A 
point imaged in any of the first four spectral bands shall be registered 
to the same point in any of the other first four spectral bands within 
.2 IFOV along-track and across-track. This registration error Includes 
the effects of all synchronizing and correcting signals developed by the 
instrument that are required to process the data into an image. Simi- 
larly Bands 5 and 7 shall be registered to Band 6 within .2 of Band 6 
IFOV, and Band 5 shall be registered to the first four bands within .3 
IFOV. These requirements shall be satisfied at least 9055 of the time 
and apply to all points along a scan line." Visual blink comparison of 
the different bands in both TM-4 and TM-5 for registration accuraoy was 
extremely good; detailed analysis was required to find subtle errors. 

Band-to-band misregistration was assessed using a correlation 
algorithm capable of estimating shifts to subpixel accuracy and which 
can be applied row-wise or column-wise. The algorithm finds the lag at 
which the correlation coefficient between any two bands is maximum by 
using cubic spline interpolation. It was found that better results are 
obtained if the registration algorithm is applied to images in which 
discontinuities are enhanced by some kind of gradient operation techni- 
que. The gradient operation used was implemented as follows; if gl cor- 
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respond to the j. uge plane of Band i, then for aoross-track measurements 
of misregistration, every row in was convolved with 

[- 110 ] 

and the magnitude of the resulting image was taken. This is equivalent 
to computing the magnitude of the gradient along the rows of g^. The 
advantage of this is that it enhanoes all the boundaries and sharp dis- 
continuities in the image, which is exactly the kind of scene features 
an observer uses when evaluating misregistration by visual inspection. 
The gradient operator was applied row-wise for estimating the registra- 
tion in the aoross-track direction or column-wise for the along-track 
direction. 

The means and standard deviations of the shifts for blocks of 
64x256 was computed using P-tape data from Illinois for both TM-4 and 
TM-5. Misregistration for the TM-4 have been previously assessed [4], 
and the possible existence of misregistration problems between the pri- 
mary focal plane Bands (Bands 1-4) and Bands 5,7 was reported. More 
recent results obtained for TM-4 by preprocessing the data with a gradi- 
ent operator Indicate that Bands 1 through 4 are registered within the 
tolerance limits with Bands 5 and 7. Bands pairs 4,1 and 4,3 gave mis- 
registration of at most .38 pixel, which is above the tolerance limit 
of .2 IFOV. This misregistration could result from the observed poor 
correlation between Bands 1,3 with Band 4 and therefore is not necessar- 


ily conclusive 
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For TM-5 we found that the first four bands , which lie in the 
primary focal plane, have a mean misregistration of, at most, .07 pixel, 
which is well within the .2 pixel specification. Results for the esti- 
mated across- and along- track registration for TM-5 are presented in 
Table 1 and Table 2, where is the mean misregistration between Bands 

i and j, and o y is the corresponding standard deviation. Overall, the 
ij 

TM-5 complies very well with all the registration requirements for both 
along- and across-track directions.. 

Striping Effects 

There are two basic types of striping effects that have been 
observed on TM-4 and TM-5 imagery. These are detector- to-deteotor 
striping and scan-to-scan striping (also called banding). The first is 
due to gain differences between the detectors in any given band, and the 
second to alternating levels of brightness related to the forward and 
reverse scanning operations. 

TM-5 detector- to-detector striping and banding was assessed by com- 
puting the mean values for 256 consecutive scan lines over a Lake Michi- 
gan scene (July 18, 1984). This area was selected because of the low 
intensity level and small variation present on the signal, facilitating 
the study of banding and striping phenomena. Figures 1 and 2 show 
results 'obtained for TM-5 Bands 4 and 7. In Bands 4 and 7, the banding 
phenomenon was clearly seen, with an average difference between scans of 
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approximately .9 and .51 binary number. In Band H the difference was 
almost one binary number, which made it easily noticeable by visual 
inspection. Striping was also noticeable in these two bands. Band 5 
showed striping formed mainly by a single detector, with a standard 
deviation of approximately .8 binary number. Also Bands 6 and 7 showed 
evidence of striping. In Band 2 an oscillation with a period of approx- 
imately 256 scan lines was observed. The differences in binary number 
observed in the TM-5 due to banding and striping are all below 1, indi- 
cating that relative calibration is very good. 

For TM- 1 !, banding and striping effects were also studied over the 
same area of Lake Michigan. In TM-*J banding was observed only on Band 
*1, while striping was present in all other bands. The striping observed 
in Bands 7,5,?, 1 was due mainly to a single detector. With the excep- 
tion of Band 7, the banding and striping deviation is below the one 
binary number range. For Band 7 striping due to a single detector was 
in some cases as big as 1.1 binary number. 

Coherent Noise Analysis 

Presence of coherent noise in the TM- 1 ! has been reported by diffe- 
rent resear cherst 1 !, 5]. The source of this noise is believed to be 
related to specific electronic or mechanical oscillations in the system. 
TM-5 data were searched for the presence of coherent noise sources. An 
average power spectrum was estimated over a block of 256x256 in Lake 
Michigan. The power spectrum measurement used is given by 
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256 

201og 10 r^ 5 DFT<U) jgiJ> 3 l<j<256 
i = l 

i 

where g^j is the image with corresponding row and column index i,j and 
w^j is a Hamming window function. Figure 3 shows the estimated power 
speotrum for TM Band 3> where the presence of three noise frequencies 
can be easily seen with corresponding periods of 12.8, 6.24, and 5.12 
pixel. The power in the TM-5 noise frequencies is nominally 5 to 7 db 
above the background noise level in Bands 2, 3, 4, 5, 7* 

Table 3 compiles the different noise periods measured over Lake 
Michigan in both the TM-5 and TM-4. The measured noise periods in the 
TM-4 and TM-5 are not the same. The observed periods were too low to be 
attributed to jitter. No in-depth attempt was made here to find the 
source of this noise. Results are presented as evidence to aid 
researchers and data processing agencies in correcting for these anoma- 
lies . 


Resolution Estimation in TM Data 

One geometric parameter of particular interest is the actual in- 
flight resolution of the sensor system which includes all environmental 
and data processing factors. The resolution is determined by the 
point-spread function (PSP) of the system. This can be estimated in 
image data by observing the measured response to scene elements of known 


shape . 
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For a space-invariant PSF, the measured image in the spatial domain 
can be expressed as a convolution of the scene with the PSF, i.e., 

g(x,y) = h(x,y)« f(x,y) 

where f(x,y) is the earth scene, h(x,y) is the overall PSF of the sensor 
system, and g(x,y) is the resulting image. 

Three scene elements that are useful for this type of analysis are; 

1. An Impulse represented by a narrow-width discontinuity 
along a row or column of the data, e.g., narrovr roads. 

2. A step function represented by an abrupt change in gray 
level along a row or column of the data. 

3. A rectangular pulse represented by a sequence of two steps 
in opposite directions along a row or column of the data. 

For our purposes, we estimated the line-spread function (LSF) in 
the direction normal to the discontinuity generated by any of the three 
suggested scene elements. If the LSF is known for many different orien- 
tation angles [2], then the two-dimensional PSF can be reconstructed. 

It was observed by experimentation with the three types of scene 
elements that the step edge was the most desirable to use in Bands 5 and 
7. In particular, agricultural field edges with high contrast fields of 
uniform texture were used. For Bands 1 through ^ the first type of 
scene element was used. To determine the LSF, subpixel sampling is 
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required sinoe the sample rate of the TM is nominally the width of the 
predicted LSF. A procedure wa3 defined to allow subpixel sampling of 
edges or roads (IEEE 1984). One critical aspect of the data extraction 
procedure is the accurate estimation of the scene element orientation, 
3ince any error here will propagate to the LSF estimate. Performance 
curves whioh relate the variance of the scene element orientation to the 
s,lgnal-to-nolse ratio were obtained, and a disoriminant funotion for 
scene element selection was developed [1 ). The subpixel scheme was 
applied to scene elements in the Iowa and Illinois data sets. A set of 
dense samples was obtained for several scene elements in all bands, 
except Band 6. For the case of field edges, the LSF was estimated in 
earlier work[6,7] using a finite sum of basis functions or by relating 
the derivative of the system response to the LSF. 

When a narrow road is used for estimating the LSF, it can be 
modeled as a narrow-width discontinuity superimposed on a constant back- 
ground, l.e., 

f <x> = Arect<^r> + B 

where t is the width of the road, A is the amplitude of the discontinu- 
ity, and B is the constant background. The resolution of the TM is 
approximately 30 meters; therefore, if rural area roads with a width of 
approximately 6 to 7 meters are used, the system response can be used to 
establish upper bounds in the system resolution. The equivalent width 
of the LSF has been previously used to characterize the resolution of 
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the TM. It is interesting to see how the equivalent width of the TM to 
a narrow road oan be related to the equivalent width of the actual LSF. 
It can be showntl] that the exaot relation is 


= « w eq <9> 



where a is going to be les3 or equal than one due to the raonotonically 
decreasing behavior of the LSF as you move away from the origin. Since 
we do not know h offhand, a is not known exactly. But, for the case of 
T = 6 meters and h given by a G&ussian function with half amplitude 
width of 1, we obtain = .995. This convinced us that narrow roads can 
be used for estimating the system resolution capabilities as long as the 
half-timplitude width of h is approximately equal to one or greater than 
one. The estimated LSF is given by a smoothed representation of the 
system response to the road. 

Scene elements for the two cases considered were extracted from the 
Iowa and Illinois scenes. Figures *1 and 5 show the raw data for Bands 3 
and 7. The resulting LSF, when cubic splines are used for smoothing, 
are shown in Figures 6 and 7- Table 4 is a compilation of the LSF 
half-amplitude widths estimates, using two types of smoothing algorithm, 
in the TM-5. The units used in this table are sampling intervals (si), 
where one si = 28. 5 meter. As we concluded in earlier work for the 
TM-*!, TM-5 resolution estimates also indicate that the actual overall 
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resolution is slightly less than the Joint electronic and optical system 
resolution of 1 IFOV. This would be due to atmospheric and data pro- 
cessing effects or to the characteristics of the scene elements used. 


Calibration of L.andsat-5 TM Thermal IR Data 

The procedures for calibration of the Landsat-4 TM thermal IR data 
and the resulting water temperature maps have been reported in a previ- 
ous paper [4]. This paper deals with the calibration of the thermal IR 

data collected by the Thematic Mapper sensor on board Landsat-5 and it 
includes a description of the procedures developed for calibration, the 
results of calibrated thermal IR data, and a comparison of these proce- 
dures and results with those obtained from the Landsat TM-4 data. 

Calibration Procedure 

The calibration of the TM-4 thermal IR data was straight forward 
because the temperatures of the Internal calibration sources (blackbod- 
ies) were known with reliability, i.e., 260 K and 320 K, which corres 

pond to the 0 and 255 digital count levels, respectively. On the other 
hand, the temperatures of the reference blackbodies for the TM-5 sensor 
were not known with. a high degree of accuracy[8]. Therefore the first 
step in the calibration procedure involved the determination of these 
two reference tempertures for the TM-5 sensor. 
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To determine these temperatures, the minimum (Rmin.) and maximum 
(Rmax.) spectral radiances stored as radiometric calibration ancillary 
records on the P-tape CCTs[9] were used in conjunction with the results 
of the integration of Planck's equation for a bandwidth defined by the 
mean system response function of all four band-6 detectors. Table 5 

p 

shows the values (in mVJatts/cm .sr. m) of the minimum and maximum speo- 
tral’. radiances corresponding to the 0 and 255 digital count levels for 
the TM-5 thermal IR band. 


Table 5. Spectral Rmin. and Rmax. for the TM-5 Thermal IR Band 


Rmin. (Spectral) = 0.1237800 inWatts/cm^.sr .um 
Rmax. (Spectral) = 1.5599560 mVJatts/cm 2 .sr .urn 


These spectral radiances were then converted to "in-band" radiances 
for a band defined by the mean system response function of all four 
detectors . The system response curves for the four detectors of both 
the Landsat 4 (protoflight ) and Landsat-5 (flight) sensors are illus- 
trated in Figure 8, and the "in-band" Rmin. and Rmax. radiances are 
given in Table 6. 
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Table 6, "In-Band" Rmin. and Rmax. for the TM-5 Thermal IR Band 


Rrain .( In-Band ) = 0.2438466 (mWatts/cm^.sr .) 

2 

Rmax .( In-Band ) = 3*0731133 (mWatts/cm .sr.) 


These Rmin. (In-band) and Rmax. (in-band) radiances correspond to the 
amount of energy emitted by a blackbody radiating at 203*2 K and 341.2 
K respectively, which should be the temperatures of the two internal 
calibration sources for the thermal IR band of the Landsat-5 Thematic 
Mapper sensor system, and they should correspond to the 0 and 255 digi- 
tal count levels of the data. 

Since the Thematic Mapper Product (P-tape) data have been radiome- 
trically corrected so that there is a linear relationship between radi- 
ance and the corresponding digital counts (DC), the following linear 
transformation was derived for the conversion of digital counts into 
in-band radiances and vice versa: 

R (In-Band) = (DC + 21.977) / 90.127 Eq. 1 

2 

Where, R = In-Band Radiance in roWatts/cm .sr. 

DC = Digital Counts between 0 and 255 
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Planck's equation was integrated between 10.00 pm and 13.00 pn at 
spectral intervals of 0.01 pm and for a range of temperatures between 
203.2 K and 3-4 1.2 K. The resulting "in-band radiances" were then multi- 
plied by the system response values (the mean for all four TM-5 detec- 
tors) given by Markham and Barker[10], and then summed over the entire 
10,00 pm and 13-00 pm temperature range, as illustrated by Equation 2. 



A + (i+1) AA 

Li 

2TTc 2 A" 5 (e hc/XkT - 1)“ 1 d A ) 

A. + iAA 
Li 


Where, 


p 

W = In-band radiance in mWatts/cm .sr. 

c = Speed of light 

h = Planck's constant 

A = Wavelength 

k = Boltzmann's constant 

T = Temperature in degrees Kelvin 

SR = System response (Markham and Barker [10]) 

A a 10 m 

AA = .01 m 

i = 1-300 


Eq. 2 


The "in-band" radiances derived from Equation 2 can be closely 
approximated by "in-band" radiances produced by a square band with lower 
and upper limits located at the 50)6 cutoff of the system response curve, 
l.e., 10. *12 pm and 12. *15 Pm, respectively. 
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As reported by Bartolucci et al.[11], the relationship between 
radianae (or digital counts) and temperature, for ranges of temperatures 
greater than 10 degrees Kelvin is not linear; therefore a quadratic and 
a cubic fit were derived to convert the TM-5 thermal IR band digital 
counts into radiant temperatures. Figure 9 shows both the quadratic and 
cubic regression curves, and for illustration purposes only, the linear 
regression curve also is included in the figure. The linear, quadratic, 
and cubic regression models are given by the following expressions 
(Equations 3» 4, and 5): 


Linear Model 

T(K) = 219.972 + 0.526 DC Eq. 3 

Quadratic Model 

T(K) = 209.831 + 0.8311 DC - 0.00133 DC 2 Eq. 4 

Cubic Model 

T(K) = 206.127 + 1.0545 DC 2 - 0.00371 DC + 6.606E-6 DC 3 Eq. 5 

Both quadratic (Eq. 4) and cubic (Eq. 5) models may be used to con- 
vert the TM-5 thermal IR band digital counts into radiant temperatures. 
However, depending upon the range of temperatures in question, one or 
the other model (either the quadratic or the cubic model) will better 
fit the actual data, as reported by Bartolucci et al[12]. 
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Temperature Mapping of Water Bodies 

To test the TM-5 thermal IR data calibration procedure, a Landsat-5 
TM scene (ID: 50139-16035) collected on July 18, 1984 over the Chicago- 
Illinois, region was utilized. First, a hierarchical classification was 
performed using data from the 1.55 m to 1.75 m middle IR band to sepa- 
rate the water pixels from all other ground cover types. Then, the 
relative thermal IR spectral response (digital counts) corresponding to 
water pixels was calibrated using Equations 3, 4 and 5* The resulting 
radiant temperatures of water were compared with reference temperatures 
measured at seven different locations (five in Lake Michigan, and two in 
the Dresden nuclear power plant cooling lake). Table 7 shows the radi- 
ant temperatures computed using the linear, quadratic and cubic models, 
the original digital counts, the reference temperatures, and the radiant 
temperatures obtained from the integration of Planck's equation without 
using regression fits. 

In Table 7 notice that the cubic model results show a better fit for 
the actual data or no-model temperatures, for a temperature range bet- 
ween 16.4° and 30 . 2 ° c. Also the radiant temperatures (no-model temper- 
ature column) and the reference temperatures (ground temperature column) 
show a discrepancy in absolute values from a minimum of 0.3° C in Lake 
Michigan, up to 2.4° C In the Dresden cooling lake. These discrepancies 
are to be expected because in this thermal IR data calibration procedure 
the atmospheric effects have not been taken into consideration. 
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Furthermore, the reference temperatures were obtained a few hours 
Defore and after the Landsat-5 TM sensor colleted the thermal XR data. 
Nevertheless, these preliminary results indicate that there is great 
potential for utilizing the relatively high spatial resolution thermal 
IR data from the Thematic Mapper sensors for measuring and mapping earth 
surface temperatures with a high degree of accuracy and reliability. 

Additional analyses of TM-5 thermal IR data in conjunction with 
more reliable measurements of reference temperatures and use of atmos- 
pheric effects correction models should be conducted before a final ver- 
dict is given for the Landsat TM thermal IR data collection capabili- 
ties. 


Summary and Conclusions 

This paper describes results of Landsats 4 and 5 sensor data qual- 
ity analysis activities which were carried out over a three-year period 
covering the Landsat 4 and 5 launches and deployment. The geometric and 
radiometric errors found in the TM-5 data were not large and their 
effect on utilization of the data is considered to be small. The per- 
formance must be considered in terms of original specifications. For 
geometry, two parameters are most significant: band-to-band registra- 
tion and large area geometric fidelity. Band-to-band registratin tests 
indicate that all bands are within or near specification. Large area 
geometric accuracy evaluation is not reported in this paper. 
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Radiometric analysis indicates that some striping and scan-direc- 
tion dependent banding still exists but with small amplitude. Banding 
is most severe in Band A with nominally a one-count difference between 
east and west scans. Several coherent noise frequencies still exist in 
the datay however, they are also of low amplitude. The total energy in 
these frequencies is very small; no known effects on data utility are 
expected. 

Resolution estimation of TM-5 detectors produced results similar to 
TM- 1 ! with IFOV values slightly larger than the specifications. This is 
expected to be cue to a combination of smoothing because of processing 
effects plus atmospheric effects. Thermal calibration activities pro- 
duced a method for obtaining a temperature for each digital count for 


TM-5. 
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Table 1. Misregistration in the ucross-traok direotion over a 256x64 


blook in Illinoia. 


Band 

1 

2 

3 

4 

5 

7 

6 

1 

0.00 

0.00 

0.00 

0.06 

0.04 

0.02 

0.32 

2 


0.00 

0.00 

0.07 

0.02 

0.01 

0.45 

3 



0.00 

0.06 

0.02 

0.01 

0.49 

4 CX ± j 3 s 




0.00 

-0.G2 

-0.02 

-0.30 

5 





0.00 

-0.01 

0.04 

7 






0.00 

-0.32 

6 







0.00 








- 


1 

2 

3 

* 'V 

5 
7 

6 


0.06 

0.07 

0.21 

0.16 

0.14 

2.76 

0.00 

0.02 

0. 19 

0.11 

0. 12 

2.71 


0.00 

0. 12 

0.08 

0.10 

2.68 



0.00 

0.12 

0.08 

1.70 




0.00 

0.07 

1.85 


0.00 1.62 
0.00 


\ 




Table 2 


TM-5 misregistration in the along-traok direction over a 64x256 
block in Illinois. 


Band 


1 

2 

3 

4 

5 

7 

6 

1 


0.00 

0.01 

0.01 

-0.06 

0.07 

0.06 

0.36 

2 



0.00 

0.01 

0.00 

0.05 

0.04 

0.67 

3 




0.00 

-0.01 

0.02 

0.01 

0.40 

4 

[ V ■ 




0.00 

0.05 

0.06 

0.14 

5 





0.00 

-0.01 

0.00 

7 







0.00 

-0.25 

6 


— 






0.00 

1 


0.00 

0.03 

0.02 

0.53 

0.17 

0.15 

1.85 

2 



0.00 

0.02 

0.23 

0.12 

0.11 

1.82 

3 




0.00 

0.07 

0.10 

0.09 

2.01 

4 

Co 3 = 

« i 




0.00 

0.19 

0.12 

1.89 

5 





0.00 

0.03 

1.79 

7 







0.00 

1.31 

6 








0.00 




Table 3. Coherent noise periods (in pixel cycle) 


using Lake Michigan data. 


Band TM- 1 ) TM-5 


1 

17, 3.12 

none 

2 

17, 3.12 

6. 24 

3 

17, 3.12 

12.8, 6. 24, 5.12 

M 

17, 3.12 

5.12 

5 

3.12 hardly noticeable 

5.12 

7 

3-12 hardly noticeable 

12.19, 5.12 

6 

none 

none 



Table 4. Results of TM resolution estimation in terms of half- 


amplitude width. 


Band 

Trig. Smoothing, VJ1/2(si) 

Spline Smoothing, W1/2(si) 

1 

1.42 si 

1.25 

si 

2 

1.20 si 

1.16 

si 

3 

1.112 si 

1.31 

si 

4 

1.50 si 

1.06 

si 

5 

1.33 si 

1.34 

si 

7 

1.37 si 

1.49 

si 


si = sampling interval 
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Figure 1. Soan-llne averaging for low-level TM-5 Band 4 data over a 

256x256 block showing scan-to-scan striping. 

Figure 2. Scan-line averaging for low-level TM-5 Band 6 data over a 

256x256 block showing detector-to-detector striping. 

Figure 3* Power spectrum of low-level TM-5 data showing noise peaks at = 
5.12, 6.24, 12.8 pixel/oycle. 

Figure Extracted road raw data from rural area in Illinois. 

Figure 5. Extracted step edge raw data from Illinois scene. 

Figure 6, Estimated LSF from step edge for TM-5 Band 7* 

Figure 7. Estimated LSF from the system response to a road for TM-5 Band 

3. 

Figure 8. System response curves for the Landsat-lJ (proto-flight) and 
Landsat-5 (flight) thermal IR band detectors (after Markham and 
Barker, 1985C10], 

Figure 9. Linear, quadratic, and cubic regression models for the rela- 
tionship between digital counts and temperature for the Landsat-5 


thermal IR band 
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Figure 6. Estimated 
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and Landsat— 5 (Flight) Thermal IR Band Detectors. 
[After Markham and Barker, 1985] 
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Piaure 9. Linear, Quadratic and Cubic Regression Models 
S for the Relationship Between Digital Counts 
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Figure 6. The Loadings or Coefficients of the Original Wavelength Bands 
on Each of the Ordered TM Principal Components of the Chicago 
0*Hare Test Site. 
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Figure 7. The Loadings or Coefficients of the Original Wavelength Bands 

on Each of the Ordered TM Principal Components of the Des Moines 
Iowa Test Site. 
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Figure 8. Percent of Total Variance Accounted for by Each of the 
Ordered TM Principal Components of the Des Moines Iowa 
and Chicago O'Hare Test Data Sets. 
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Figure 11. 


Image of Band 5 (1.55 - 1.75uin) of Landsat 
IV TM of the Dresden Power Plant Test Site. 



Figure 12. Image of Band 6 (10.4 - 12.5yra) of Landsat 
IV TM of the Dresden Power Plant Test Site. 
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Figure 34. Temperature Map of the Dresden Nuclear Power Plant 
Test Site Derived from the Landsat-IV TM Data. 
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Abstract 

Land sat-4 thematic mapper and multispectral scanner data were analyzed 
to obtain information on data quality and information content. Geometric 
evaluations were performed to test band-to-band registration accuracy. The- 
matic mapper overall system resolution was evaluated using scene objects 
which demonstrated sharp high contrast edge responses. Radiometric evalua- 
tion Included detector relative calibration, effects of resampling, and 
coherent noise effects. Information content evaluation was carried out 
using clustering, principal components, transformed divergence separability 
measure, and numerous supervised classifiers on data from Iowa and Illi- 
nois. A detailed spectral class analysis (multispectral classification) was 
carried out on data from the Des Moines, Iowa area to compare the informa- 
tion content of the MSS and TM for a large number of scene classes. 

Introduction 

Landsat-4 thematic mapper and multispectral scanner daua were analyzed 
in this study to obtain Information on data quality ard information con- 
tent. Geometric evaluations were performed to test band-to-band registra- 
tion accuracy. Thematic mapper overall system resolution was evaluated 


* This work is sponsored by the National Aeronautics and Space Administra- 
tion under NASA Contract NAS5-26859. 


using 3cono objects which demonstrated oharp high oontraat edge responses. 
Radiometric evaluation included detector relative calibration, coherent 
noiae analysis and the effects of resampling. Information content evalua- 
tion waa carried out uaing clustering, principal components, and the trans- 
formed divergence separability measure on data from Iowa and Chicago, Illi- 
nois. A detailed spectral class analysis (multispcctral classification) 
wa3 carried out on data from the Des Moines, Iowa area to compare the 
Information content of the MSS and TM for a large number of scene classes. 
This work is being carried out as part of a Landsat-4 data quality evalua- 
tion study. 

T est Data Sets 

The aata sets used in this study were gathered over the Des Moines, IA 
area on Sept. 3, 1982 and the Chicago, IL area on October 25, 1982. Both A 
tapes which are radlometrlcally corrected and P tapes which are both 
radlometrically and geometrically corrected were acquired for these scenes. 
The scene numbers were 40049-16264 for Iowa and 40101-16025 for Illinois. 
Preliminary test scenes over Detroit and Arkansas were also available. 

Comparison of A and P Tape Data 

To determine the effect of the geometric correction process (resam- 
pling and gray-level interpolation) on the radiometry of the new resampled 
pixels when converting the data from an A- tape (30m x 30m pixels) to a 
P-tape (28.5m x 28.5m pixels) format, the means and standard deviations for 
homogeneous and heterogeneous areas on the ground were calculated using \ he 
data from both the A and P tapes. A qualitative comparison of both sets of 
statistics Indicated that no change (radiometric degradation) was caused by 
the non-zero order (cubic convolution) Interpolation process. 
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Furthermore, to determine whether the cubic convolution interpolation 
had affected the structure of the data in feature space, a set of data from 
the A and P tapc3 for the same area on the ground was clustered into 16 
cluster classes. Then the resulting 16 cluster classes from the A and 
P-tape data sets were merged together into a single 32-class file. Subse- 
quently, the pair-wise spectral separability for the 32 classes was com- 
puted uring a Transformed Divergence algorithm!! ]. The minimum transformed 
diverge^ je between cluster classes from the TM A tape and corresponding 
second cluster classes from the TM P tape i3 22 and the maximum is 262. 

Since the transformed divergence (D ) measure rarges from zero to 2, 000 — 

T 

where a value of 2,000 indicates that the two classes in question are com- 
pletely separable (different) and any D T value between zero and 500 indi- 
cates that the pair of classes is not separable (very similar) — it may be 
concluded that the 16 cluster classes obtained from the A and P tapes are 
essentially equivalent. Thus, these results corroborate the previously 
stated conclusion that the radiometry of the resampled pixels from the P 
tape has not been significantly affected by the geometric correction pro- 
cess. Consequently, they meet GSFC specifications for the Landsat-D sys- 
tem, i.e., resampling shall not introduce radiometric degradation of more 
than 1 quantum level[2]. 

Registration Evaluation 

The first TM scene, which was of the Detroit area, contained several 
geometric and radiometric errors which did not appear in subsequent scenes. 
The overall visual quality of TM data is extremely good; detailed analysis 
is required to find subtle errors. Band-to-band registration was evaluated 
using a correlation algorithm capable of estimating shifts to subpixel 
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acouraoy, The algorithm Implements the correlation coefficient and user 
cubic interpolation to estimate the peak to subpixel accuracy. The means 
and standard deviations of shifts for 32 sets of lines in TM Imagery of the 
Iowa and Arkansas sites were computed. Results indicate that the primary 
focal plane bands (1, 2, 3, 4) are registered to .1 pixel, which is well 
within the .2 pixel specification. Bands 5 and 7 are also well registered 
to within the .2 pixel specification. 

Problems were discovered, however, in the registration of the first 
band set (1, 2, 3, 4) to Bands 5 and 7. Approximately .6 pixel misregis- 
tration was observed between these bands, which exceeds the specifications 
of .3. The thermal IR band could not be numerically correlated due to the 
great difference in the nature of the image relative to any of the other 
bands, so a visual blink comparison was used on a digital display unit. A 
four 30 M pixel line misregistration was observed in both Arkansas and 3 owe 
thermal data. 

The MSS band-to-band registration was similarly evaluated and shifts 
were found to be less than the .2 pixel tolerance for the A data, but 
slightly greater than the tolerance for P data. The maximum shift of .27 
was between Bands 1 and 2. 

Detector-to-Detector Calibration 

Thematic mapper detector-to-detector calibration was evaluated by com- 
puting mean values of sets of scan lines spaced 1 6 lines apart In A-tape 
data. Results are presented in Table 1. Differences as large as 1.3 digi- 
tal count were observed, but the typical deviations from the overall moan 
were less than .5 digital count. These results indicate that striping 
should be minimal in TM data. This i3 verified by visual inspection of the 


imagery. 
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Similar tests were run on MSS data. The results are presented in 
Table 2 in terms of maximum deviation from overall means and variances. 
The values are all below one bin, indicating that relative calibration is 
very good. 

Coherent Noise Analysis 

Visual inspection of MSS imagery over Lake Michigan revealed signifi- 
cant noi3e patterns in all bands. These wavelike patterns had a period of 
to 15 samples with higher frequency patterns superimposed. The magni- 
tude of the noise is 4 to 5 counts peak-to-peak and it exists in all bands. 
It was first assumed that this noise was a jitter effect. Fourier analysis 
was performed on samples of the noise patterns and several significant fre- 
quencies were observed. The largest had wavelengths of 3.31, 4.09 and 15.4 
pixels which agreed with the visual inspection. The frequencies were all 
too high, however, to be caused by jitter. A sampling, power supply, or 
data processing is suspected as the cause. 

No visual noise effects were apparent in the thematic mapper data. 
However, Fourier analysis of low level Lake Michigan data revealed two 
coherent noise frequencies in some of the bands. Figure 1 contains a plot 
of the power spectrum of data from band 4 detector 1 and two noise frequen- 
cies are apparent. One 13 at a wavelength of 17 pixels and the other is at 
3.12 pixels. The power in these frequencies is 10 to 20 decibels (10 to 10C 
times) above the background noise level, but relative to the power in the 
bright areas of the scene the noise power is very low. The amplitude of the 
3.12 pixel noise signal is 24 decibels and the low frequency power in the 
Sept. 23, 1982 Iowa TM data is in the 60 decibel range, thus the noise is 
more than 30 decibels less than the image power. This corresponds to a 
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standard deviation of nominally 30 times less than the signal. Band H stan- 
dard deviation i3 21 counts for typical sites in the Iowa data, thus the 
amplitude of the noise is nominally 0.6 count. These noise frequencies 
exist in all detectors of band 4. This noise wa3 observed in bands 1,2,3 
also, but not in bands 5,6 and 7. 

The MSS and TM noist frequencies are nominally the same suggesting a 
common cause, but no explanation is advanced here. These results are pre- 
sented as evidence to aid researchers and data processing agencies in cor- 
rection and refinement of the sensor and data processing stream to aid in 
assuring maximum quality of Landsat-^ data. 

Resolution Estimation in TM Data 

One geometric parameter of particular interest is the actual in-flight 
resolution of the sensor system which includes all environmental and data 
processing factors. Th-* resolution is determined by the point-spread func- 
tion (PSF) of the system. This can be estimated in image data by observing 
the measured response to scene elements of known shape. 

The measured data can be expressed in the spatial domain as a convolu- 
tion of the scene *ith an overall point-spread function: 

g(x,y) = h(x,y) * f(x,y) 
where f(x,y) is the earth scene 

h(x,y) is the overall point- 
spread function of the 
sensor system 

g(x,y) is the resulting Image 
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Given g(x,y), we wish to determine h(x,y). To do this, some deterministic 
element of the input f(x,y) must be known or assumed. Although the theory 
can take into account the two-dimensional nature of the element, the ini- 
tial experiments have been limited to the one-dimensional case. If the 
overall PSF is separable, l.e., if h<x,y) can be written a3 a product 
h<x)h(y), then this approach provides a direct estimate of the two compo- 
nents. Otherwise It generates cross sections through the two-dimensional 
PSF along the x and y axes. 

Three scene elements that would be useful for this type of analysis 

are: 

1. An impulse represented by a narrow-width discontinuity along a row or 
column of the data . 

2. A step function represented by an abrupt change in gray level along a 
row or column of the data. 

3. A rectangular pulse represented by a sequence of two steps in opposite 
directions along a row or column of the data. 

It was observed by experimentation with the three types of scene ele- 
ments that the step edge was the most desirable to use. In particular, 
agricultural field edges for high contrast fields with uniform texture were 
used. To determine the PSF, subpixel sampling is required since the sample 
rate of the sensor Is nominally the width of the PSF. A procedure was 
defined to allow subpixel sampling of edges. 

Because of the orbital inclination of Landsat and the propensity of 
man to arrange linear features, such as roads and field boundaries, in the 
cardinal compass directions, It is generally found that there is a spatial 
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displacement In the scene coordinates of the linear elements from one row 
or column to the next. This ha 3 the desirable effect of providing a fine 
grid of samples of the system response when values from adjacent rows or 
columns are combined after correction for the spatial shift of the scene 
element. The procedure for combining the data is quite straightf orward and 
can be illustrated as follows for a north-south road. 

The coordinates of the peaks in the row data corresponding to the road 
are determined for a sequence of N rows. These data are then fitted with a 
least-squares straightline , providing an analytical expression for the road 
coordinates. The x-coordinates in each row are then modified by subtract- 
ing from them the least-squares estimate of the road location in that row. 
This converts the data to a coordinate system in which zero is the road 
center. Because of the small angular difference between the sensor coordi- 
nate system and the road direction, the change in road coordinates from row 
to row is only a fraction of the pixel spacing and so represents a sampled 
response to the scene element corresponding to a subpixel translation. By 
combining the data from a number of rows, a set of finely sampled data is 
found. These data can be graduated using splines or other smoothing func- 
tions to give an average response function from which to estimate the sys- 
tem point-spread function. 

The subpixel sampling scheme was applied to field edges in the Iowa 
data set. A dense set of samples was obtained for several edges in TM 
Bands and 5. The PSF was estimated in earlier work[3] using a finite sum 
of basis functions. The method used here is based on the following obser- 


vations s 
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The 3 cene element can be modeled as a step superimposed on a constant 
background, l.e., 


f (x ) = A + B u(x) 

where A is the background and B is the amplitude change across the boun- 

0 

dary. The output is then given by 

g(x) = [A + B u(x)] * h<x) 

Taking the derivative of both sides of this equation gives 
g'(x) = (A+B u (x) ) * h(x) 

= B5 (x) * h(x) 

= B h(x) 
h(x) = g'(x) 

The PSF is therefore the derivative of the measured step function. The 
measured subpixel sampled step data were smoothed to produce a well-behaved 
derivative. If this had not been done, the derivative would have been 
excessively noisy. Two methods of smoothing were used: (1) Trigonometric 

polynomials, and (2) cubic spline functions [4]. 

A simulation experiment was conducted using an ideal edge with Gaus- 
sian noise to test the smoothing methods. A test PSF was convolved with 
the noisy edge and the result was smoothed and the PSF :d.dth was determined 
by differentiation. It was observed that the spline smoothing tended to 
underestimate the width by 10^ and the trigonometric smoothing overesti- 
mated by 15%. The results were used as correction factors when the method 
was applied to real data. 
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Field-edge data for two cases were extracted from the Iowa 3oene and 
the derivative PSF estimation scheme was applied. The resulting estimated 
point spread function for Band 4 is shown in Figure 2. The results for the 
two types of smoothing are listed in Table 3. The width measure is the 
half amplitude of the estimated PSF. These results indicate that the 
actual overall resolution is less than the optical system resolution. This 
would be due to atmospheric and data processing effects or to the charac- 
teristics of the field edge. The effect of smoothing may also be signifi- 
cant but the simulation experiment provided a correction which is believed 
to compensate for this. These results are basfed on a limited sample under 
one atmospheric state and further estimation research must be carried out 
to determine the variability of the resolution estimates. 

Scan Angle Effects 

The effect of scan-look angle on the recorded digital response value 
was studied using the TM A Tape from Iowa (Scene ID: 40049-1 626*1 ) . Average 
digital values calculated from 3 blocks each of dimensions 16 lines (detec- 
tor) by 10 columns were obtained for both the forward and reverse scans; 
each set of 3 blocks was selected at intervals of 100 columns for each of 
the reflective bands. Figure 3 shows a plot of the average digital response 
values for the forward and reverse scans for Band 1. This figure shows a 
general decline in average response across-track, with a total drop of 
about 4 to 5 digital counts. Some of the other channels displayed similar, 
although somewhat lower and less distinct drops in digital counts across- 
track. Much of the variance across-track was due to scene variation and 
both linear and multiple regressions failed to produce significant equa- 
tions for all bands (Table 4), but a general decline in response does seem 
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to exist for some if not all of the bands, as indicated by a negative slope 
obtained for all regressions from both reverse and forward scans for all 
bands. No significant difference between the forward and reverse scans for 
any band was found for the Iowa TM A Tape using an alpha level = 0.05. 

Dimensionality of TM vs. MSS 

The spectral dimensionality of the TM sensor is of great interest 
since it potentially represents nearly double the information over the MSS. 
The dimensionality was investigated by applying the principal components 
transformation to test data. 

One of the major advantages of using orthogonal linear transforma- 
tions, such as the Principal Component or Karhunen-Loeve transformation, Is 
their ability to compress overall data variance in a multifeature space 
onto a relatively few transformed orthogonal axe3. This essentially sepa- 
rates non-random variance (information content) from random variance 
(noise) while, concurrently, uncorrelating the transformed axes in such a 
way that any information redundancy (due to inter band correlation) is eli- 
minated. Such transformations therefore allow a compression of the data 
into a fewer number of dimensions while retaining a maximum amount of sig- 
nificant Information content and removing much of the random variance or 
noise from the data. 

Principal components were generated for both the MSS and TM P-tape 
data sets from the Chicago O'Hare test site. Statistics used in calculat- 
ing the principal components were generated from data samples of the origi- 
nal MSS’ and TM data sets from every fifth line and fifth column and all 
bands. Since the sum of the eigenvalues Is equal to the trace of the ori- 
ginal covariance matrix, l.e., the total variance, the importance or per- 
cent of total variance explained by eigenvector (X^) is given by: 
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tr S 


where * eigenvalue i or characteristic root 
of the 1th component 

tr S = trace of covariance matrix S or 
total variance. 

Therefore, each of the eigenvalues of the ordered components divided by the 
sum of the eigenvalues represents the amount of total data variance or 
Information content accounted for by each eigenvector. Tables 5 and 6 list 
the eigenvalues and the corresponding amount of data variance that is 
accounted for by their respective eigenvectors for both the MSS and the TM 
data sets. As shown in these tables, the first two principal components of 

the MSS P-tape account for almost 97% of the total data variance. This is 

also evident from the images of the ordered MSS principal components (Fig. 
iJ). The first two components contain a significant amount of scene con- 
trast and data structure, while the last two MSS components contain mostly 
random noise. 

In the TM data set, the first two ordered components together account 

for 90$ of the data variance, the first three account for 97. 1 J5$i and the 

first four components account for almost 99$ of the data variance. The 
images of the TM principal components also show that the first four compo- 
nents contain a significant amount of scene contrast with a sequential 
decrease in the lower-ordered components (Fig. 5). 

From both the listed eigenvalues of the MSS and TM data sets, it is 
possible to conclude that the first two components of the MSS and the first 
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four components of the TM data sets contain a maximum amount of significant 
Information. It could be possible , therefore, to conclude that for this 
data set of the Chicago O'Hare area, the TM data contain more significant 
information, i.e., a higher significant dimensionality, than the corres- 
ponding MSS data set. The significant dimensionality as defined by a prin- 
cipal component analysis may be substantially different for other data sets 
and other ground cover types. However, these results do imply that for the 
same area, the TM may contain more significant information than the corres- 
ponding MSS. 


Pr.' nolpal Component Comparison; 

(Chicago O'Hare and Des Moines Iowa TM Data) 

In order to investigate both the significant dimensionality and the 
Importance of the individual TM bands for different data sets, a principal 
component analysis was performed on the Des Moines Iowa TM data (Scone 
ID: 40049-16264) and compared with the principal components analysis of the 
Chicago O'Hare TM data set discussed above. 

Since the coefficients of a particular principal component correspond 
to the direction cosines of the original bands to the new component [ 5 ], it 
is possible to describe the relative direction of the new component axis in 
the feature space with respect to the original bands. The greater the 
coefficient of a particular band to a specific component, the larger the 
cosine and thus the smaller the angle subtended between the original fea- 
ture axis (band) and the new component axis, i.e., the more the new compo- 
nent axis lies in the direction of the original feature axis (band). In 
this way It Is possible to describe the relative influence or "pull , ' of the 
original bands on each of the new component axes. 
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Using the coefficients of the higher-ordered principal components, 
l.e., those components which incorporate most of the total data variance, 
it is possible to describe which of the original bands contain most of the 
significant variance or information content for a particular data set. Fig- 
ures 6 and 7 graphically show the loadings or coefficients of the '’hicago 
O'Hare airport and the Des Moines, IA TM data sets, respectively. In addi- 
tion, Figure B graphically shows the percent of total data variance 
accounted for by each eigenvector (principal component) of the Chicago 
O'Hare and Des Moines, IA TM data sets. 

It is interesting to compare these two data sets in that they demons- 
trate the importance of different bands for different data sets of varying 
cover type. The first eigenvector or component of the O'Hare data set shows 
relatively equal weighting of all of the original reflective bands 
(Fig. 6). By comparison, the greater amount of vegetation (primarily agri- 
cultural crops and forest areas) associated with the Des Moines TM data set 
has resulted in a significant amount of the total variance (accounted for 
by PCI, Fig. 8) of the data to be contained on Band H and, to a lesser 
extent, Band 5. The second component of the De3 Moines 7“ data set contains 
more information from the rest of the bands. It is Interesting to note that 
component four of both the Chicago O'Hare and Des Moines TM data sets has 
the greatest contribution from band 6 (thermal IR) (Fig. 6 and, 7), which is 
clearly evident from a comparison of the two images, i.e. component H and 
the original band 6 (Fig. 5). This may imply that although the thermal IR 
band may not contain a significant amount of total data variance (the 
fourth component of both data sets describe only 1 —2% of the total data 
variance), the Information or variance may be distinctly unique from the 
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rest of the bands. In addition, although most of the significant variance 
was incorporated on different bands for these two data sets (as described 
by the first eigenvector for each set), both sets indicate a significant 
dimensionality of between. 3 and iJ dimensions (Fig. 8). 

Spectral Class Analysis in TM and MSS Data 

A detailed spectral analysis was conducted of thematic mapper and MSS 
data for an area near Des Moines, Iowa from the September 3, 1982 data set. 
Data were utilized from 7 blocks distributed throughout the area which 
included agricultural, forest, suburban, urban, and water scene types. The 
blocks were processed using a clustering algorithm to produce up to 18 
cluster groupings for each block. Each cluster class was then identified 
with a ground-cover class using aerial photography and maps of the area. 
The cli!f k ,r classes from each of the 7 blocks were inspected with regard to 
separabililty, means, variances, and number of points in each class were 
either deleted or pooled with spectrally similar clusters[1 93 . The separ- 
ability measure used in the transformed divergence function or proces- 
sor [6] measures the statistical distance between classes based on class 
means and covariance matrices . The measure has a maximum value of 2, 000 
and the minimum of 0. Spectrally, very close classes will typically have 
values as low as 50 to 500. 

For the TM data, initially 9*1 classes were defined and the pooling and 
deleting process reduced these to *J2 final spectrally separable classes. 
Table 7 lists these classes. The MSS data were then analyzed using the 
same clustering and merging sequence. The number of separable classes in 
the MSS is 21, half of the TM result. This result is considered to be a 
very significant indicator of the dimensionality of TM relative to MSS. 



PAGE If. 


The MSS olaao occurrences are indicated in Table 7. The maximum divergence 
values of any one class with respect to all others also were much less fo* 
the MSS classes relative to TM. Table 8 contains the maximum and average 
transformed divergence values for the *J2 spectral olasses and for the best 
subsets of TM spectral bands. It should be noted that the best speotral 
band for any combination of Bands 1 through 7 is the first middle IR band 
( 1.55-1. 75ym) . The next best band is the near IR (0.76-0. 90pm), followed 
by the red band and then the thermal IR. The best combination of four 
bands includes one from each of the four regions of the spectrum {visible, 
near IR, middle IR, and thermal IR). Table 9 contains the minimum and 
average transformed divergence values for the best combination of MSS 
bands. 

The high average divergence indicates that the 21 spectral classes 
found in the MSS were about as separable as the TM classes. However, there 
was twice the number of equally separable TM classes. This is considered 
to be the most significant result of the spectral analysis. 

A final test using a preliminary set of test data was carried out 
using the small amount of ground truth available. A set of 5,615 TM and 
1,376 MSS pixels containing forest, corn, soybean, soil, water, and urban 
classes was extracted from the TM and MSS data where the cover classes were 
known or could be inferred from aerial photography. The overall correct 
recognition was 95.7/6 for the TM classification using all 7 bands, 92. 6 % 
for the TM using a subset of the best four bands, and 67 - for the MSS 
classification using all H bands and a per-point, Gaussian Maximum Likeli- 
hood (GML) classifier. The results are listed In Table 10 along with the 
amount of CPU time required to classify each data set. 
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An additional classification using all 7 TM bands wa3 performed on the 
test data using a contextual algorithm, SECHO (Supervised Extraction and 
Classification of Homogeneous Objects), in order to demonstrate the effec- 
tiveness of such contextual algorithms over per-point algorithms for use 
with the 30 meter resolution TM data. The contextual or per-fleld algorithm 
SECHO first divides the scene to be classified into homogeneous fields and 
then classifies these fields using an extension of the GML algorithm [73* 
SECHO Incorporates the fact that since cover classes are more likely to 
occur in homogeneous areas larger than one pixel in size (i.e. larger than 
30 meters), adjacent pixels are highly correlated, with the degree of cor- 
relation diminishing with an increasing distance between the pixels [73. 
Thus SECHO assigns an analyst-specified threshold value, below which adja- 
cent pixels will be grouped into a homogeneous field. Statistics for these 
fields are calculated and compared to the original cover class statistics 
and a "homogeneous field" is classified as a unit into that class which it 
most closely resembles. 

These "test" fields are limited in number of pixels and so are not 
really evaluating how representative the final spectral classes are of the 
entire scene but rather how separable the classes are. Deletion of certain 
spectral classes (e.g., Corn2), due to low separability, resulted in much 
confusion of corn with trees in the MSS but not the TM for both the best 
four and for all seven TM bands. Also, the resolution of TM actually 
allowed "purer" cluster classes to be defined since smaller areas (e.g., 
beaches, roads) were distinct. These results along with the listed GPU 
time for classification Indicate that a subset of the best four TM bands 
incorporates the advantages of the higher spatial resolution (e.g., "purer" 
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cluster cla33ca) over the equivalent MSS data without simultaneously incur- 
ring as substantial an increase in computer time required for classifica- 
tion as with all seven bands. In addition, these test results support the 
use of contextual algorithms, such as SECHO, over per-point algorithms for 
use with TM data. An earlier study using simulated TM data [8] demonstrated 
that higher classification performances could be expected with the SECHO 
classifier for the 30 meter TM data, especially in those areas where the 
cover class field sizes are relatively large, relative to the scanner FOV, 
e.g., most agricultural areas. 

Temperature Mapping of a Cooling Pond and 
Thermal Plume from a Nuclear Power Plant 

The use of thermal Infrared (XR) data in conjunction with reflective 
multispectral data, obtained from aircraft altitudes, has been proven to be 
an effective means for increasing the classification accuracy of earth sur- 
face features [9,10]. 

Also it has been demonstrated that radiometrioally calibrated thermal 
IR data obtained from airborne platforms can be used to produce accurate 
temperature maps of water bodies [11,12,13,14]. 

Accurate temperature maps of water bodies have been obtained also from 
calibrated thermal XR data gathered in 1973 by the SKYLAB SI,-2 SI 92 scanner 
system[15]. However, for almost ten years the remote sensing data users 
have not had available high spatial resolution thermal XR data (in contrast 
to the coarse spatial resolution data gathered by meteorological satel- 
lites) acquired from spaceborne platforms, until the TM thermal data were 


recently collected. 
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In order to determine the radiometric quality of the TM thermal data 
for temperature mapping of surface water, a test site was selected within 
the area covered by the TM scene (Scene ID: *10101-16025) gathered over 
Illinois. This site was. chosen because it includes a surface water body 
with a large range of temperatures, i.e. a cooling pond for the Dresden 
nuclear power plant and the Junction of two rivers: the Kankakee and the 
Des Plaines Rivers which once merged form the Illinois River, 

The Dresden power generating station is located approximately 50 miles 
southwest of Chicago near Morris, Illinois. The station houses three 
nuclear reactors which together have a maximum nameplate generating capac- 
ity of 1,656,630 Kw. The water which is taken from the Kankakee River and 
is used to cool the stream condensers is channeled into a man-made lake 
which is designed tu dissipate heat. This 1300 acre lake has a dike in its 
center which directs the water flow up, around, and down the lake in an 
eight- mile loop that takes two-and-a-half days to complete. During its 
course the waier is cooled by natural evaporation and is either recycled 
back through the station or discharged into the Illinois River. A schematic 
diagram of the power plant, channels, cooling lake, and the three rivers is 
shown in Figure 9. 

Because the radiant temperature of the various cover types in the test 
site is a function not only of the kinetic temperature of the materials, 
but also a function of the Intrinsic emissive properties (emissivity) of 
the objects, and to avoid differential emissivity complications, it was 
decided to perform the calibration of the TM thermal Id data corresponding 
to only those pixels representing water, which has an emissivity close to 
that of a perfect radiator or blackbody. 
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To accomplish the thermal calibration of water only, a heirarchical 
classification of the TM data had to be performed, l.e., one of the middle 
IB (1.55 - 1.75 wm) bands was utilized to discriminate the water pixels 
from all the other cover types using the LARSYS layered classifier 
[12,16,173. Once all the water pixels were separated from everything else, 
their thermal IR responses (relative digital counts) were converted to 
radiant temperatures. Figure 10 shows the decision tree utilized to carry 
out the heirarchical or layered classification and calibration of the water 
bodies present in the test site. Figures 11 and 12 show the images of the 
test site corresponding to the 1.55-1.75ym and 10. t|-12. 5ym bands, respec- 
tively. 

The conversion of the relative thermal IR responses into radiant temp- 
erature measurements was accomplished using a non linear calibration func- 
tion derived specifically for the TM thermal IR band and the range of temp- 
eratures of the two internal calibration reference blackbodies. 

The non-linear calibration function was derived through the Integra- 
tion of Planck's equation for a spectral band covering the 10. 4-12. 5pm 
spectral range, and for a range of temperatures between 260 K and 320 K 
[18] at Increments of 0.1 K. The resulting in-band radiances (in 
Watts/cm 2 -sr ) were then plotted versus the corresponding temperatures, as 
illustrated in Figure 13. Note in this figure that if one assumes a linear 
relationship between the emitted energy from a blackbody and its tempera- 
ture for a range of temperatures in excess of 10 K (the range between Ihe 
cold and hot internal calibration sources of the TM system is 60 K), errors 
on the order of 4 K (4°C) would be introduced. 

The ron-linear calibration function is given by Equation 1 below: 


PAGE 21 


T(°C) = -12.5809 + 0.2917D - n. ’'00233 D' <D 

where: T = temperature in degrees Celsius 

D = relative digital counts for the TM thermal IR 
band 

If a linear interpolation was performed, the linear calibration func- 
tion would be represented by Equation 2: 

T(°C) = -13 + 0.23529D (2) 

The pixels of water at the intake (Kankakee River) have a relative 
thermal IR response equal to 95 digital counts, and those corresponding to 
the water at the entrance of the cooling pond equal to 135 digital counts. 

If these two relative thermal responses (digital counts) are substituted in 
Equation 1, one obtains the following radiant temperatures: 

Digital Count Radiant Temperature 


95 

13. 1 °C 

135 

22. 6° C 


Since at the time the TM data were recorded (October 25, 1982 at 16:02 
GMT) reference (ground truth) data were not being collected, It was not 
possible to assess reliably the accuracy of the radiant temperatures der- 
ived from the TM thermal IR data. However, the personnel from the Dresden 
nuclear power plant keep a record of the temperature of the water at the 
entrance to the cooling lake. The temperature recorded for that date was 
71°F or 21.7°C, which is a difference of approximately 1°C with respect to 
the radiant temperature derived from the TM thermal IR data. Figure 1H 
shows a temperature map of the Dresden power plant test site. 
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Although neither the temperature of the water recorded at the power 
plant nor the precise temperature of the Internal calibration reference 
blackbodies in the TM system was reliable, the results obtained during this 
preliminary assessment of the TM radiometric data quality are very encour- 
aging. 

Summary and Conclusions 

Landsat-4 thematic mapper and multispectral scanner data quality was 
evaluated using data from several frames over midwestern United States. 
Four categories of evaluation were presented: 

1. Basic radiometric and geometric characteristics. 

2. Overall TM system resolution. 

3. Spectral dimensionality. 

4. Spectral class separability. 

5. Use of TM thermal IR data for temperature mapping of water bodies. 

The overall quality of the TM data is very good with only a few error 
conditions, which in the case of band-to-band misregistration can be cor- 
rected by ground processing. 

Detector- to- detector calibration was very good, such that striping 
should be minimal in TM. MSS demonstrated more noise problems, primarily 
due to hig£i frequency wave patterns. In general, the MSS appeared somewhat 
noisier than previous Landsat MSS' data. 

The resolution estimation experiments resulted in relatively wide 
half-amplitude widths for the point-spread function of the TM. The reason 
for this is not understood; it is likely due to a combination of system and 
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atmospheric effects plus the uncertainty in the actual shape of the step 
edges being used as a test input. 

Dimensionality was investigated using principal components analysis. 
The results indicate two . significant dimensions in the MSS and up to four 
in the TM. The spectral analysis was the most revealing by producing twice 
a3 many separable classes in TM than in MSS. The result was due both to 
the greater number of spectral bands and to greater spatial resolution of 
the TM. More classes, such as beach, roads, field edges, and other small 
features, were separated out in TM where they are mixtures in MSS and in 
many cases not spectrally distinct. These results further verify the high 
quality of the TM data and promise significant increase in usefulness of 
these data in earth resources applications. 

The radiometric quality of the TM thermal IR band was studied using 
both a linear calibration function and a non linear calibration function 
derived from the Integration of Plank's equation for a spectral band cover- 
ing the 10. - 12.5vm spectral range, and for a range of temperatures bet- 
ween 260 K and 320 K at increments of 0.1 K. These preliminary results sup- 
port the use of the non linear calibration function with the thermal IR 
band for effective use of the thermal band for temperature mapping of water 
bodies. 
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i*ble ' 

DlANNCL-to-CHAHHCL CALIBRATION t VALUAT10S 
Deviation of Detector Mein fron Orind he*n lor VLBSTCR CO,, IOWA *A“ DATA 

5e»n Line Sequence 


Dp- 

tector I 

2 

3 

4 

5 

6 

7 

B 

9 

10 

11 

12 

13 

14 

15 

16 

Bind 
















1 -.08 

-.4 

.14 

-.04 

.29 

.32 

-.38 

.44 

.35 

.07 

-.1 

-.1 

.49 

-.24 

.2 

-.29 

2 .S3 

-.41 

.35 

-.47 

.73 

-.08 

.06 

.16 

.2 

.03 

-.41 

-.55 

-.55 

-.02 

.42 

.05 

3 .45 

-.13 

.19 

-.05 

.1 

-.77 

-.14 

-.2 

-.23 

-.11 

.21 

-.14 

-.39 

.04 

.43 

.76 

4 .2 

.4 

.42 

.18 

1.0 

.9 

.64 

-.1 

-.77 

-.94 

-.65 

-.59 

-.38 

-.22 

.08 

-.17 

5 .12 

-.15 

.11 

-.48 

.05 

-.02 

.49 

-.25 

-.05 

.04 

-.23 

-.44 

-.44 

-.*4 

1.29 

.33 

6 -.29 

0.0 

.32 

.04 

.69 

.25 

.24 

.36 

-.21 

.61 

.17 

.02 

-.39 

.18 

.26 

-.35 
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Table 2. 


Maximum Deviations of Mean and Variance 


of any Detector for MSS Bands > 


WEBSTER CO. 


CHICAGO 


A Data 


A Data 


Band 

Mean 

Std Dev. 

Mean 

Std. Dev 

1 

-.28 

.43 

.2 

.16 

2 

.30 

co 

oj 

• 

.2 

.25 

3 

.82 

.29 

.13 

.23 

H 

.4H 

.2 

.11 

-.12 
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Table 3- 

Results of TM Resolution Estimation 
Using Step Differentiation. 
Average of 3 Test Regions 
"A" Tape Data 30M Pixel Spacing 
TRIGONOMETRIC SMOOTHING SPLINE SMOOTHING 




W (Plx.) 

Corr. 

W (M) 

Knots 

W (Pix.) 

Corr. 

W (M) 

Band 

4 

1.60 

1.44 

43.2 

7 

1.13 

1.31 

39.3 

Band 

5 

1.62 

1.46 

43.8 

7 

1.29 

1.49 

44.7 


